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Abstract 

This paper presents a calculation of the quasilocal energy of a generic FRW model of 
the universe. The results have the correct behavior in the small-sphere limit and vanish 
for the empty Milne universe. Higher order corrections are found when comparing 
these results to classical calculations of cosmological energy. This case is different from 
others in the literature chiefly in that it involves a non-stationary spacetime. This fact 
can be used to differentiate between the various formulations of quasilocal energy. In 
particular, the formulation due to Brown and York is compared to that of Epp. Only 
one of these is seen to have the correct classical limit. 

1 Introduction 

In general relativity, as other areas of physics, the notion of energy is of both conceptual 
and practical importance. Unlike other areas of physics however, a sensible and universal 
definition of energy has been difficult to come by. Difficulties in defining energy in the 
context of general relativity became apparent from the very beginning. Einstein's pseudo- 
tensor formulation, as well as the many variations that have appeared since [TJ[2], are not 
entirely satisfactory because they are not generally covariant. The global definitions of 
energy, such as the Bondi or ADM formulations, also seem inadequate because they are 
not localized and must necessarily refer to the entire spacetime|^ A reasonable compromise 
seems to have been reached with the notion of 'quasilocal energy' - energy that is defined for 
an extended but finite region of spacetime. There is a plethora of definitions of quasilocal 
energy, many of which are thoroughly reviewed in [2] . 

The virtues and weaknesses of these definitions can be revealed by application to specific 
spacetimes. The Schwarzschild and Kerr spacetimes are common testing grounds for the 
various notions of quasilocal energy. However the fact that these spacetimes are stationary 
limits their discriminating power. Non-stationary examples can discern among the various 
notions of quasilocal energy in a way that stationary examples cannot. The primary purpose 
of this paper is to examine the case of a generic Friedmann-Lemaitre-Robertson- Walker 
universe, commonly referred to as FRW cosmology in much of the literature. Aside from 
its practical applicability, this non-stationary example reveals in a new way how various 
notions of quasilocal energy differ from one another. Due to limitations of time and space, I 
have chosen to focus on only two formulations of quasilocal energy among the many possible 
ones. The first is the formulation of Brown and York [4 and the second is the closely related 
formulation of Epp [5]. I will comment on other formulations only in passing, although they 
certainly deserve greater attention. 

From a theoretical point of view, the Brown- York formulation is appealing because it 
is based on the rigorous Hamilton-Jacobi theory of classical mechanics. Although one may 

^See however [3] for a mass function that 'localizes' the ADM energy, albeit in a computationally in- 
tractable manner. 



question the applicability of this theory to general relativity, the ends ultimately justify its 
use. So this formulation will be my starting point in Sec. [2j Calculation of the Brown- York 
quasilocal energy, Eby, for the the FRW model is carried out in Sec. |2.1| where it is seen 
to be contradictory to physical expectations. In Sec. [3] the Epp quasilocal energy, Ee, is 
introduced, and the calculation for the FRW model is repeated in Sec. |3.1[ It will be shown 
that in the appropriate limits, this second result provides a reasonable and expected value 
of energy for an FRW universe. 

2 Quasilocal Energy a la Brown and York 

To begin with, the notation will be as follows. M denotes the spacetime manifold. I 
assume that M has topology M x S, where S is a fixed three-dimensional manifold. This 
assumption allows the foliation of M into spacelike leaves, Sj, each of which represents space 
at one instant in time. To be more precise, I assume that there exists a diffeomorphism 
/ : E X E ^ M. Then for t e M and a f {t, a) represents one point of the leaf E^. Note 
that / = (/°, f^,P, P) , and that it can be used to define a vector field = df^/dt on M. 
This vector field specifies the direction of time evolution at each point. R denotes a region 
in St for some fixed i, and denotes the boundary of this region, = dR. B denotes 
the manifold that results when Q. evolves in time according to the vector field t^^. In other 
words, B denotes the manifold that results by moving points of along the integral curves 
of V^. Note that B will have the topology M x fij^ See Figure [l] for a simplified depiction of 
the geometry. 




As in 13], each of the manifolds described above is endowed with a metric: g^i, is the 
metric on M, while /ly , 7ij and Cab are the induced metrics on E, B and O respectively. Note 
the use of indices, where ^i,v = 0,1,2,3 are used to index coordinates on M, i,j = 1,2,3 
are used (somewhat confusingly) to index coordinates on either S or _B, and a,b= 1,2 are 
used to index coordinates on f2. The unit normal to S in M is denoted by m^, and the 
extrinsic curvature with respect to this normal is denoted by K^^. The unit normal to B in 
M is denoted by rt^, and the extrinsic curvature with respect to this normal is denoted by 
The unit normal to f2 in S is denoted by n^, and the extrinsic curvature with respect 
to this normal is denoted by fc^. Table [l] summarizes some of this notation. 

The main result of |4] is that a generalized energy-momentum tensor, r*-', may be de- 
fined which describes the total energy, i.e. energy of matter content as well as the gravi- 
tational energy, inside an arbitrary region. The motivation for this definition comes from 
the Hamilton- J acobi formalism in classical mechanics. If Sdas denotes the action calculated 

■^The notation here closely resembles that of Brown and York 4_, with the exception that where I use B 
they use ^B, and where I use O they use B. Epp [5J uses similar notation, with the exception that where I 
use Q he uses S. 



2 



Manifold 


Dims. 


Metric 


Normal 


Ext.Curv. 


M 


4 








S 


3 








B 


3 


7ij 






n 


2 




rii 


h 



Table 1: Summary of Notation 



for a classical solution, i.e. a solution satisfying the equations of motion, then momentum 
can be expressed as the derivative of S^ias with respect to coordinates on the boundary. 
Similarly, Brown and York define the generalized energy-momentum tensor r'-' as: 

ii _ 2 SSclas 



Here Sdas is the total action for gravity and matter calculated for a classical solution, and 
the variation is with respect to the metric on the boundary B. Calculating the functional 
derivative in ([T|) shows that 

" 2 „ 



where tt*^ = dS/Sjij is the canonical momentum conjugate to the boundary metric. The 
generalized energy-momentum tensor is used to define the surface energy density e, surface 
momentum density j", and surface stress density as follows: 

e = UiUjT^-' (3) 
f ^ -atu.T^^ (4) 
s =a.^(TjT-'. (5) 

With the appropriate choice of coordinates, the generalized energy-momentum tensor can 
be written as: 



r 



J 

The quasilocal energy inside the region bounded by f2 is then defined as: 



(6) 



E= / d^x^/\c7\ e. (7) 

To simplify equation ([t]) note first that, using the same mathematics as in the ADM 
formalism, one can express the canonical momentum as: 



|7| (67'^' - e'^") , (8) 



where Qij is the extrinsic curvature of i? as a submanifold of M. The surface energy density 
can then be written as: 

2 1 

e = u.Ujt'^ = =MjUj7r'^ = ^^"^^j (©7'^ - ©'^ ■ (9) 

Decomposing Qij into components tangent and normal to S, it can be shown that e = 
k/SirG, where k is the trace of the extrinsic curvature of as a submanifold of Sj^ This in 

^It should be emphasized that k = a'^^kab where kab is the extrinsic curvature of Q viewed as a subman- 
ifold of E. One can define another extrinsic curvature for f2 viewed as a submanifold of B. This point will 
prove to be important. 
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turn allows one to write equation (|7| as: 



E ■ 



inG 



\a\ k. 



(10) 



As with most definitions of energy, a reference point must be chosen. Since it is generally 
believed that the energy in Minkowski spacetime must be zero, the standard approach is to 



use this spacetime as the reference. Equation ( 10 1 is thus modified to yield the Brown- York 
quasilocal energy: 
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(11) 



where quantities with an overbar refer to the reference spacetime and those without refer 
to the physical spacetime. Specifically, f2 is a suitably chosen 2-manifold in Minkowski 
spacetime, and k is the trace of the extrinsic curvature of that 2-manifold. There remains 
some ambiguity in the choice of fl. In [J] it is suggested that fl ought to be an isometric 
embedding of f2 into Minkowski spacetime. This prescription can be applied to a variety of 
simple examples with very good results. For example, it can be applied to static, spherically 
symmetric spacetimes |4], and the results are in agreement with expectations. However, 
when applied to the non-stationary case of a Robertson- Walker spacetime, the results do 
not seem to be correct as demonstrated in the following subsection. 



2.1 FRW Cosmology a la Brown and York 

FRW cosmology is described in comoving coordinates {t, r, 6, (j>) by the Robert son- Walker 
metric: 



1 — Kr^ 



r'^dO'^ + sin^ ( 



(12) 



where k g {—1, 0, +1} is proportional to the spatial curvature and a(t) is the dimensionless 
scale factor. 

To simplify the calculation, region i? is chosen to be spherically symmetric. The bound- 
ary of this region, fi, is then a sphere of comoving radius r at time t. Adopting the natural 
constant-t foliation of M, one may calculate the extrinsic curvature of in S: 



diag 







-r-\/ 1 — nr^ a{t) , — r 1 — nr'^ sin^ 6 a{t) 



(13) 



The first term in (|ll|) is then calculated to be 

sin 0a^{t) I - 



E, 



phys 



1 

8^ 



2vr 



r a(t) 



-^jy/l- Kr^a{t). (14) 



Alternatively, one can take a more direct approach and calculate the generalized energy- 
momentum tensor: 



1 



inG 



= diag 



\/l — r^/T^^~Kr^ a{t) r\/l — Kr^ sin^ 0a(i) 



(15) 



AnGrait) SttG SttG 

The t-t component of this tensor is the surface energy density e which, when integrated over 



fi, provides the same result as (14) 



Ephys ^ [ e = --^rv^r^-KT^ a(t). 



(16) 
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To calculate the reference term, I choose Minkowski spacetime as the reference spacetime 
and denote all quantities in the reference spacetime with a bar over them. So Minkowski 
spacetime, M, is coordinatized in spherical coordinates by = {t, r, 9, (p) , and has the 
metric 

g^, = diag[-l , 1 , f2 , f^sin^^ . (17) 

Before calculating the reference term, one must find an isometric embedding, ip : M, 
of the boundary CI in Minkowski spacetime. I use the following ansatz which preserves 
spherical symmetry: 

i;{e,(t>) = {tit,r) , r(t,r) ,0,0). (18) 

The functions t(t, r) and f(t, r) are to be determined by the requirement that the metric be 
preserved, i.e. that the embedding be isometric. This embedding defines Cl. The metric on 
Cl is the induced metric, i.e. the pull-back of Minkowski metric to the submanifold Cl: 

a = tp*{g) ^ (Tab = diag [f^ , sin^ 9] (19) 

To enforce the requirement of isometry, the metrics on and Q are equated, aab = (^ab, with 
the result that f{t,r) = a{t)r. Note that equating the metrics does not determine t{t,r), 
although this ambiguity will not affect any of the following results]^ Now the second term 



in ( 11 ) can be calculated: 



1 /• ...T _2 . .f-2\ 1 . 1 



Subtracting ( 20 1 from ( 14 1 , the Brown- York quasilocal energy of a Robertson- Walker 



spacetime is found to be: 

Eby 



(l - y/l - Kr2) a(t). (21) 



The obvious deficiency in this result is that for the case of a flat FRW model, k = 0, the 
quasilocal energy is zero regardless of the matter content of the model. Clearly something 
is missing here. 

Quasilocal energy is proportional to the extrinsic curvature. However, for an arbitrary 
2-manifold embedded in a 4-manifold, the extrinsic curvature has two components. In other 
words, being of codimension two, has two linearly independent normal vectors each of 
which has an extrinsic curvature associated with it. A natural question then is to ask which 
of the two components should be used to calculate quasilocal energy. Above, is viewed as a 
submanifold of S, and its normal in S is used to calculate the quasilocal energy. However, 
one could have just as well viewed f2 as a submanifold of B and used its normal Ui in B to 
perform the same calculation but with different results. To put it differently, the calculation 
above depends in a critical way on the choice of foliation. A different foliation of M would 
mix the two components of extrinsic curvature, giving an entirely different expression for the 
quasilocal energy. The static spacetimes considered in ^ have a natural foliation for which 
the second component of extrinsic curvature vanishes and thus makes no contribution to 
the quasilocal energy. For a Robertson- Walker spacetime, this is no longer the case. What 
is missing from the above calculation is a proper treatment of both components of extrinsic 
curvature. 



^If one were interested, there is a way to determine t{t, r) as well. Instead of embedding Q in M to 
obtain f2, one could choose to embed B in M, and then take the appropriate cross-section to obtain Q. 
This procedure would provide an explicit expression for t(t,r). Calculation of the quasilocal energy however 
would be the same in either case. 
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3 Quasilocal Energy a la Epp 



The basic idea behind the notion of quasilocal energy, at least in the Brown- York and 
Epp formulations, is that the energy contained in a closed surface is related to the extrinsic 
curvature of that surface. Therefore to begin with, a precise definition of extrinsic curvature 
is needed. 

Let M denote an arbitrary manifold of dimension m, and TpM denote the tangent vector 
space of M at a point p. Let TpM be spanned by linearly independent vectors {e^}, where 
H— 1, - ■ ■ ,m. Let g denote the metric and V the covariant derivative on M. 

Let S denote a submanifold of M of dimension s, and let TpS denote the tangent vector 
space of S at the point p. Let TpS be spanned by linearly independent vectors {^i}, where 
i — 1, • • • ,s. Let NpS denote the normal vector space of S at the point p. Let NpS be 
spanned by linearly independent vectors {ua}, where a = I, - ■ ■ ,m — s. 

An arbitrary vector v g TpM can be decomposed into two components, one normal and 
the other tangent to S: 



w - ^ g{v,na)na 



(22) 



normal 



tangent 



A similar thing can be done for the covariant derivative. Given two arbitrary vectors u,v G 
TpM, the covariant derivative can be decomposed into normal and tangent components: 



(23) 



normal 



tangent 



This decomposition can be written as Vut; = K(u,v) + VuU where K is the second funda- 
mental form and V is the induced covariant derivative on S. The second fundamental form 
is in fact defined by this decomposition [6]: 



K {u,v) = ^5 (V„t;, Ua) Ua- 



(24) 



The induced covariant derivative is the projection of V„w onto S, and is the unique covariant 
derivative that is torsion-free and compatible with the induced metric on 5 [6] . 

The mean extrinsic curvature vector H is defined as the trace of the second fundamental 
form: 



(25) 



where the vectors £ TpS are tangent to the submanifold S. The mean extrinsic curvature 
scalar is defined as the magnitude: 



\H\ = VH -H 



(26) 



To connect with the standard physics literature, it might be helpful to express some of 
these results in 'indexed' notation. The second fundamental form, mean extrinsic curvature 
vector, and mean extrinsic curvature scalar can be expressed respectively as: 



a=l 

s 



\H\ = ^g^.mn- . 



(27) 

(28) 
(29) 
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A minus sign appears in ( 27 1 when the derivative operator is moved to a neighboring term, 
and a projection operator should be used as appropriate when restricting these quantities 
to the submanifold S. 

Note that according to its definition, the second fundamental form is a tensor of type 
(1,2). In the special case where S is of codimension one, the second fundamental form is 
defined simply as K (u,v) = g (V„w, n). In this case, since there is only one normal vector, 
it can be omitted from the definition, making the second fundamental form a tensor of type 
(0,2). In this context, the second fundamental form is also called the extrinsic curvature. 

It can be shown [S] that K = (— i£„(7)", where £ denotes the Lie derivative. It can also 
be shown for u,v TpS and V compatible with g that K{u,v) = — g(w,V„n). This is the 
basis for writing, as one often sees in the literature, Kij — —ViJij. Also in this context, the 
mean extrinsic curvature vector is really just a scalar and reduces to the usual definition of 
mean extrinsic curvature found in the study of surfaces embedded in M'^ for example. 

Adapting the above notation to the case at hand, M becomes the four-dimensional 
spacetime, and S becomes the two-dimensional boundary surface fi. Being of dimension 
two and codimension two, has two tangent vectors {^i,^2}i ^'^^ ^wo normal vectors 
{n'^,n2}- The Epp definition of quasilocal energy can be expressed in terms of the mean 
extrinsic curvature scalar. The energy inside fi, without the reference term, is: 

Evhys = -^ jj^x^\\H\. (30) 

A reference point must also be defined. As before, Minkowski spacetime is used as the 
reference, and the Epp quasilocal energy becomes: 

EE^E^nys'Eref^-^ jj^x^\ \H\ + ^ J_^d' X VW\ \H\ . (31) 

Epp's derivation of this result was motivated by an analogy with the formula E'^ = 

from special relativity, and so proceeded along different lines than that presented here. 

However as Epp notes 5 , the final result can be expressed in terms of the means extrinsic 



curvature scalar exactly as in (31 1 



It is instructive to compare the Epp energy (31 1 to the Brown- York energy (111. The 
main distinction of the Epp energy is that it is independent of the foliation. This is so 
because in | both components of the extrinsic curvature are accounted for and combined 
in a way that is invariant under a change of basis for the normal space NpS. Since the 
choice of foliation is tantamount to the choice of observers, this amounts to saying that the 
Epp energy is invariant under boosts of observers. For this reason, it is tempting to refer 
to Epp's definition as quasilocal mass rather than energy. However, other properties of this 
quantity compel Epp to refer to it as an 'invariant quasilocal energy' rather than mass [5j. 

For stationary spacetimes, the Brown- York energy calculated for stationary observers is 
identical to the Epp energy. To put it differently, suppose that for a stationary spacetime one 
chooses the natural constant-t foliation, where t is the coordinate along the timelike Killing 
vector. Then, the Brown- York energy calculated for this foliation is equal to the Epp energy. 
For example, both formulations results in the same expression for the quasilocal energy of 
Schwarzschild spacetime: 

In the limit r ^ oo, this result reduces to M which is exactly the ADM energy. 

The real distinction between these formulations manifests itself when considering non- 
stationary spacetimes where the natural choice of foliation can result in non-vanishing ex- 
trinsic curvature in both components. In the next section, this is demonstrated explicitly 
in the context of FRW cosmology. 
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It should be mentioned that the Epp quasilocal energy (31) is similar to the definition 
ofTered by Wang and Yau in j^- The Wang-Yau construction requires a foliation and the 
choice of an admissible time function r : S ^ K. (See [7] for precise definitions.) Although 
their approach is considerably more complicated, it can be shown that with the natural 
foliation and a constant time function, t : S ^ tg, their approach and that of Epp give the 
same result when applied to the FRW universe. 

3.1 FRW Cosmology a la Epp 

In this subsection I again calculate the quasilocal energy of an FRW universe, but this time 
using the Epp definition (31). As in section 2.1 consider the Robertson- Walker metric and 
choose to be a sphere of comoving radius r at time t. The normal space of O is spanned 
by two unit vectors which, in coordinates (t,r,9,(f)), are: 



< = (-1,0,0,0) 



(0: 



ait) 



,0,0). 



(33) 



These vectors are used to calculate the second fundamental form. The tangent space of 
is spanned by the following two unit vectors: 



0,0, 



1 



a{t)r 



,0 



0,0,0, 



1 



a{t)r sin t 



The mean curvature vector and its magnitude are: 

' -2a{t) -2(1- 



\H\ = 



\ a{t) ' ra2(i) 

2v/l-Kr2-r2 a?{t) 
r a{t) ■ 



, 



The first term in (31 1 is then calculated to be: 



E.. 



phys 



(34) 

(35) 
(36) 

(37) 



To calculate the reference term in (31 1, the same procedure as in section 2.1 can be used. 
Specifically, Minkowski spacetime is again the reference spacetime, and Cl is defined by the 
same isometric embedding (18). The normal space of Cl is spanned by two unit vectors 
which, in coordinates (t, r, 9, (f)) i ^rs: 

= (-1,0,0,0) , = (0,1,0,0). 

The tangent space of f2 is spanned by the following two unit vectors: 



1 



0,0,^,0 
r 



^2 = 0,0,0, 



1 



r sm ( 



The mean curvature vector and its magnitude are: 

-2 



iJ^ = ( , 
\H\ = 



, 



The second term in (31 1 is then calculated to be 



ref 



<rxJ\a\ \H\ 



G 



G 



ra{t). 



(38) 
(39) 

(40) 
(41) 

(42) 
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Finally, subtracting ( 42 1 from ( 37 1 , the Epp quasilocal energy of a Robertson- Walker space- 
time is found to be: 



phys 



— Eref — 



G 



nr^ — r'- a 



(43) 



This new result, as opposed to (21 1, does not generally vanish for a flat, k = 0, FRW 
model. The quasilocal energy depends on both the scale factor and the rate of expansion. 
These quantities are determined by the Friedmann equations and depend on details, such 
as matter density and equation of state, of the FRW model. 

Some obvious features of this result can be checked immediately. Consider the case of a 
flat, K = 0, static universe, a{t) = 1. In this case the FRW model reduces to the Minkowski 
spacetime, for which the quasilocal energy (43 1 is zero, as expected. Of course, equations 
(43 1 and (21 1 are equivalent for static universes. The difference becomes apparent when the 
scale factor is allowed to evolve. So consider the case where k = —1 and a{t) ~ t, which 
describes the Milne universe. It is well known that the Milne universe is just a portion of 



Minkowski spacetime expressed in different coordinates |8]. Using (43) the quasilocal energy 
of the Milne universe is zero, as expected for any portion of Minkowski spacetime. Equation 



(21 1, on the other hand, would not have given the correct result. 



3.2 Comparison to Newtonian Cosmology 



A further test of the reliability of the result ( 43 1 is to compare it with a classical calculation of 



energy in Newtonian cosmology. It is expected that in the case of a flat, matter-dominated 
universe, and in the small-sphere limit, the quasilocal energy (43) should reduce to the 
Newtonian results. 

For a flat universe (k = 0) dominated by (non-relativistic) matter, the Friedmann equa- 
tions [H] imply that a{t) — y^2Aa~''^{t), where A = ^^p^a^. Here p„ and are the matter 
density and scale factor at some initial time. Using this, along with k = 0, in equation (43 1 
results in: 



^E^-^r a{t) 



1 - Vl-2Ar2 a-i(t) 



(44) 



It is expected that in the small-sphere limit, i.e. as r — > 0, relativistic cosmology should 
resemble a Newtonian universe. To examine this limit, (44) should be expanded about 
r = 0: 



Ee = ^a{t) 



Aa-\ty + ^A^a~^{ty + O(r^) 



47r 



+ ^-^plata-'{ty + 0{r^). 



(45) 



It is encouraging to see that the first term in this expansion corresponds to the rest mass 
energy of a spherical ball of gas. To correctly compare this result with the Newtonian 
result, it should be noted that in FRW cosmology r denotes a comoving coordinate, while 
in Newtonian cosmology r usually denotes a fixed coordinate. As a consequence, in FRW 
cosmology the proper distance from origin is cr = a(t)r, while in Newtonian cosmology the 
proper distance is simply a — r. Equation (45), expressed in terms of the proper distance, 
is: 



Ep 



47r 



(46) 



Equation (46) can now be compared to the Newtonian result. A brief description of 



Newtonian cosmology can be found in [lOi |9j , while more thorough accounts can be found 
in [TTJ [T5J [13]. The appendix in this paper provides the necessary background for the 
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purpose here. Equation (64 1 from the appendix gives the total classical energy, including 
kinetic, potential and rest mass energy, inside a sphere of radius r: 



2n 



(47) 



where fc is a constant proportional to the mechanical energy, and a„ (t) is the Newtonian 
scale factorj^ In terms of the proper distance, this becomes: 



E , — 



An 



^0 NO N ^ ' 



27r 



^0 NO N ^ ' 



(48) 



The difference between the quasilocal energy (46 1 and classical energy (48 1 is then: 

Ai? = Ee — E^ias 
An 



po [ala-%t) - al^a-^t)] 



(49) 



The constant k in Newtonian cosmology plays a role akin to k in FRW cosmology. The 
case k — corresponds to an expanding gas in which every particle travels at exactly its 
escape velocity, allowing the gas to expand asymptotically to infinity. To compare with a 
flat FRW universe I choose fc = for which the Newtonian scale factor a„ (t) and the FRW 



scale factor a(t) coincide exactly. In this case, the difference (49 1 reduces to: 



AE 



8n^G 



plala-'^{t)a^ + 0{a^). 



(50) 



Thus this calculation of quasilocal energy suggests that there are corrections to the classical 
energy appearing at order O(cr^). One can be even more explicit by substituting for the 
scale factor which, for a flat matter-dominated universe, is a{t) ^ t"^^^. The correction (50 1 



then becomes AE ^ t + 0{a'^). This then suggests that although there are corrections 
to the classical energy, they tend to vanish quickly with time. 



4 Conclusion 

Defining a universal notion of energy has been a challenge since the inception of general 
relativity. Today a myriad of definitions exist. Although it is generally agreed that energy 
ought to be a quasilocal quantity, it is not at all clear which of the many quasilocal definitions 
is the correct one. Some even question whether a correct definition exists. Perhaps, like 
so many important questions in physics, this one ought to be settled by experiment or 
observation. 

The field of cosmology has proven to be a fertile testing ground for many of the ideas 
in general relativity. So it seems appropriate to look to cosmology for clues to the nature 
of energy. The primary goal of this paper has been to calculate the quasilocal energy of 
an FRW universe, the so-called standard model of cosmology. For this purpose two specific 
definitions of quasilocal energy were chosen. It has been shown that in FRW cosmology, or 
more generally for non-stationary spacetimes, a clear distinction arises between the Brown- 
York and Epp formulations. It seems that Epp's invariant quasilocal energy (or mass) is 
more appropriate for describing what we normally think of as cosmological energy. For 
non-stationary spacetimes, there is little precedence, and therefore the validity of equation 



(43 1 can be examined only in special limits. The comparison with Newtonian cosmology 



shows that the results here are reliable at least to lowest order. The presence of higher order 



distinguish between the FRW scale factor, a(t), and the Newtonian one, a^(t), because these two do 
not necessarily coincide. 
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corrections to Newtonian cosmology is of course expected. It would be interesting to examine 
how these corrections manifest themselves in cosmological data. It is conceivable that the 
presence of this additional energy may affect structure formation in the early universe, alter 
the rate of expansion, or manifest itself as dark energy [TJ . 

This work can be extended in several directions. The most immediate extension would 
be to consider other formulations of quasilocal energy. The constructions by Hawking |15j , 
Kijowski [TB], and by Liu and Yau [T7j are similar in structure to the two cases studies here. 
Several special cases of the FRW model, such as de Sitter (dS) and anti-de Sitter (AdS) 
spacetimes, deserve closer examination. In fact, using appropriate coordinates, the metrics 
for these spacetimes can be written in the form of the Robertson- Walker metric. So the 
result (43 1 is already valid for dS and AdS so long as the coordinates are carefully considered. 
The thermodynamics of these spacetimes is dependent on energy, and thus provides another 
interesting application of quasilocal energy. 



Acknowledgements 

I would like to thank Steven Carlip for offering advice and guidance without which this 
paper would not have been possible. This work was supported in part by the Department 
of Energy under grant DE-FG02-91ER40674. 



Appendix: Newtonian Cosmology 



Here I provide a review of Newtonian cosmology sufficient for deriving equation ( 47 1 . 



The setting for Newtonian cosmology is flat, 3-dimensional, euclidean space. I will use 
spherical coordinates denoted by r = (r, First, it is assumed that space is filled 

with dust of homogeneous density p{t) out to some large radius R{t) . Note that both the 
density and the radius depend on time, indicating that the ball of dust may be expanding or 
contracting. Second, conservation of mass is assumed, and in particular that the continuity 
equation holds: V • J = —dp/dt. Third, it is assumed that, inside the radius R{t), Hubble 's 
law holds: v {r,t) = H{t)r, where H{t) is the Hubble parameter. 

Hubble's law can be expressed as r = Hr. This equation can be solved to obtain: 



"(t) = roexp 



H [t') dt' 



to 



(51) 



where Tq = r (^o) for some initial time to. Using the definition J = pv and the homogeneity 
of density, the continuity equation can be expressed as pV • v = —p. Using Hubble's law in 
this last equation, it can be modified to 3pH = —p. which then can be solved to obtain: 



p{t) = po exp 



-3 I H{t')dt' 

J to 



(52) 



where po = p(to)- The Newtonian scale factor a„(t) is defined by H{t) — d„(<)/a„(i)Q 
This equation can be solved to obtain: 



fljv (*) = a«o exp 



H {t') dt' 



(53) 



where 



(to). This last equation can be used to replace H(t) with a^^it) in equations 



"Here r is a fixed coordinate, not to be confused with the comoving coordinate used in FRW cosmology. 
'^The scale factor in Newtonian cosmology is not necessarily the same as in FRW cosmology, and is thus 
distinguish here by the subscript N. 
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(511 and (52): 



r{t) 



(54) 



Pit) 



Po 



(55) 



Consider a test mass m located at a distance r < R from the origin. The force on this 
test mass, according to Newtonian mechanics, is the gravitational force of the mass enclosed 
inside the radius r; the mass outside that radius has no effect on the test mass. So the force 
on test mass m is: 

F=—^, (56) 

where M = ^nr^p{t) is the total mass inside radius r. Writing this equations as r = ~^^^rp, 
and then using equations (54 1 and (55 1 to eliminate r{t) and p(t) in favor of aj^{t), results 
in: 

««W = -^PoaLa:'(i). (57) 



Equation (57 1 can be integrated to obtain 



8ttG 



' JVO N ^ ^ 



(58) 



where k = ^^p„a^ — is a constant. This equation can be rearranged and written in 

6 "0 NO NO ^ O 

a form that resembles one of Friedmann's equations: 

2 



SttG 



6 a 



(59) 



To discuss energy in the context of Newtonian cosmology, consider again the test mass 
m located at distance r < R. From equation ( 56 1 , the gravitational potential energy of mass 
m is calculated to be: 



GMm 



AttG 



-pr m. 



(60) 



The kinetic energy of the same test mass is: 

rrn ^ 2 ^ Tt2 2 

i„, = -mv ~ -mil r . 
2 2 

Defining the mechanical energy of mass m as Um = Tm + Kn results in: 



TT 1 W2 2 47rG 2 

C/™ = -m.H r - 



pr m. 



(61) 



(62) 



Assuming uniform density, the total mechanical energy inside a sphere of radius r < R can 

Idn^G 



be found by integrating (62 1 over the volume of the sphere: 

AttG 



U = 



pH^r^ 



15 



-p'^rK 



(63) 



To find the total energy, one should add to mechanical energy the rest mass energy as well. 
So the total classical energy, including kinetic, potential and rest mass, is: 

167r2G 



An 

y 

47r 



2n 



pr 



2^5 



-pH^r- 



15 



-p^r' 



2tt 



— p„a^ a ^r^ H p„a^ a ^r^ 



IGtt^G 
15 



p^a^ a-^r^ 

f^Q NQ N 



4tt 



Po^o^n'^' 



2tt 



PoKo^^Nr^k, 



(64) 
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where in the second line I have used equation ( 55 1 , and in the third Hnc I have used the 
definition of k. 

It should be pointed out that k is proportional to the mechanical energy of a unit mass: 



fc = -2(^^^ U. (65) 

The mechanical energy is a conserved quantity in Newtonian cosmology, which explains 
why k is referred to as a constant. Most treatments of Newtonian cosmology [TTJ [T3j |9] 
consider only the case where U = and thus fc = 0. This corresponds to an expanding 
gas in which every particle travels at exactly its escape velocity, allowing the gas to expand 
asymptotically to infinity. The case fc < corresponds to an expanding gas with speeds less 
than the escape velocity, resulting ultimately in contraction of the gas under the infiuence 
of gravity. The case fc > corresponds to an expanding gas with speeds greater than the 
escape velocity, resulting in eternal expansion. 
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